Reentrant Condensation of DNA induced by Multivalent Counterions 
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A theory of condensation and resolubilization of a dilute DNA solution with growing concentration 
of multivalent cations, is suggested. It is based on a new theory of screening of a macroion by 
multivalent cations, which shows that due to strong cation correlations at the surface of DNA the 
net charge of DNA changes sign at some small concentration of cations A'o- DNA condensation 
takes place in the vicinity of A'o, where absolute value of the DNA net charge is small and the 
correlation induced short range attraction dominates the Coulomb repulsion. At N > No positive 
DNA should move in the oppisite direction in an electrophoresis experiment. From comparison 
of our theory with experimental values of condensation and resolubilization thresholds for DNA 
solution containing Spe*"*", we obtain that A^o = 3.2 mM and that the energy of DNA condensation 
per nucleotide is 0.07 ksT. 



PACS numbers; 77.84. Jd, 61.20.Qg, 61.25Hq 



I. INTRODUCTION 

In the last several years there has been a revival of 
interest in the phenomenon of DNA condensation with 
multivalent cations. The reason for this interest is the 
general effort of the scientific community to develop ef- 
fective ways of gene delivery for the rapidly growing field 
of genetic therapy. The DNA compaction should be fast, 
effective, easily reversible and should not damage the 
DNA double helix. All of these conditions are fulfilled 
in DNA condensation with multivalent cations, such as 
CoHex'^+, naturally occurring polyamines Spd'^+, Spe''"'" 
and their analogs which are known to bind to DNA in 
the predominantly nonspecific electrostatic manner. The 
DNA condensates obtained this way are indeed closely 
packed arrays of parallel DNA strands. It was shown 
that the helical structure of the B-DNA is not perturbed 
within such condensate, and that the reaction is easily re- 
versed by the addition of monovalent salt, or simply dilu- 
tion of the solution with water. Cations with larger and 
more compact charge are more effective in condensing 
the DNA. This also suggests an electrostatic mechanism 
of the DNA condensation with multivalent cations. Dur- 
ing about twenty years of research, a significant amount 
of information on DNA condensation has been accumu- 
lated. For long DNA, as the concentration of Z-valent 
cations grows condensation happens abruptly at some 
critical concentration, Nc which depends on the charge 
of cations and the concentration of monovalent salt n. A 
comprehensive review of the experimental and theoretical 
results for Nc can be found in Rcfs. 0J^. 

The intensive study of the last few years revealed com- 
pletely new features of DNA condensation with multi- 
valent cationg3i3. It was discovered that when the con- 
centration of cations grows far beyond Nc to some new 
critical value N^ ^ Nc, DNA dissolves and returns to 
the solution. This reentrance condensation behaviour is 
schematically shown on Fig. 1. For a long DNA and 
small n both transitions are very sharp and the ratio of 



Nd/Nc can be as large as 10''. Remarkably, the decon- 
densation threshold Nd is almost totally independent on 
the monovalent salt concentration, n. On the other hand, 
the condensation threshold, Nc, grows with increasing n. 
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FIG. 1. Schematical illustration of the reentrant condensa- 
tion. The fraction x of DNA molecules in solution is plotted 
as a function of logarithm of the cation concentration 

It has been understood for some time that, due to cor- 
relatioiis between multivalent cations at the DNA sur- 
faceLTtJ, two DNA molecules experience a short range 
attraction, which can lead to condensation. Monovalent 
ions are much less correlated and do not provide any at- 
traction. Conventional explanation for the condensation 
threshold Nc is that it is just the bulk concentration of 
Z-valent cations, N = Nj., at which they replace monova- 
lent ones and begin to produce a short range attraction. 
However, from this point of view it is very difficult to un- 
derstand why at large N = Nd, DNA molecules go back 
to the solution. 

The above explanation also does not take into ac- 
count the net charge of DNA—jlEhe nonlinear Poisson- 
Boltzmann equation predictsEflo a value for the net 
linear charge density of DNA, 77*, which includes the 
bare charge of DNA and the charge of cations bound 
at the very surface of DNA with energy larger than ksT. 
This charge is negative and does not depend on N. At 



1 



N > Nr, two situations are possible. The energy of the 
Coulomb repulsion can be smaller than the energy of the 
short range attraction. Then DNA condenses at N — Nr 
but never dissolves back. If the energy of the Coulomb re- 
pulsion is larger than the energy of short range attraction, 
condensation does not happen at all. Both possibilities 
contradict experiments and therefore reentrant conden- 
sation can not be explained in the Poisson-Boltzmann 
aproximation. 

In this paper we propose an explanation for the reen- 
trant condensation based on a new theory of screening of 
macroions by multivalent cations, which emphasizes the 
strong correlations of multivalent cations at the surface of 
DNA. We explain why the condensed phase exists only 
within the limited range, Nc < N < Nd and calculate 
both Nc and Nd- 

It was shown in Ref. |l^,|l^, that the strong repulsion 
between multivalent cations leads to their strong lateral 
correlations. The resulting strongly correlated liquid of 
multivalent cations at the surface of DNA has a large 
negative chemical potential, which describes the addi- 
tional purely electrostatic binding of cations to the sur- 
face. This, in turn, leads to an exponentially small con- 
centration of cations, A^O: above the surface. (Aq > A^, 
if n is not unrealistically large). According to Ref. 17 
when the bulk concentrations of multivalent cations, A^ 
grows above A^^, the net negative charge 77* decreases in 
absolute value and crosses zero at A^ = Aq- At A^ > Aqi 
the net charge becomes positive and continues to grow 
with A^. This effect is called charge inversion. It is worth 
noting here, that the chargeJnversion is not a result of 
some specific chemosorbtionE2l, but rather a direct conse- 
quence of purely electrostatic interactions of multivalent 
cations. 

At A^ = Aq there is no Coulomb repulsion at all, so 
that the short range attraction dominates and leads to 
condensation. This is actually the optimal situation for 
condensation. It is obvious that there is a range of A^ 
around Aq where attraction still dominates and DNA 
condences. Condensation thresholds Nc and Nd then are 
determined by the condition that the energy of the short 
range attraction is equal to the energy of Coulomb repul- 
sion of negative and positive DNA molecules respectively. 
This means that the concentration A'o is located inside 
the window Nc < N < Nd (see Fig. 1) and the width of 
the window grows with the strength of the short range 
attraction. 

Below we present an analytical theory, which calculates 
both critical concentrations Nc and Nd in terms of the 
two main physical parameters of the system: boundary 
concentration of multivalent cations Aq and the energy, 
£, of the DNA binding within the DNA bundle per nu- 
cleotide. 

At a small concentration of monovalent salt, n, the con- 
densation threshold Nc calculated in our theory is larger 
than the replacement concentration A^^- This means that 
the replacement of monovalent cations happens before 
the condensation and by itself can not lead to the con- 



densation while Coulomb forces are still strong. On the 
other hand, if n is so large that Nr is larger than the 
value of Nc calculated in our theory (which assumes that 
replacement has happened), the condensation actually 
happens only at A^ = A^r- Thus, a large concentration of 
monovalent ions, n, acts as a "curtain". It eliminates a 
part of the window predicted by our theory. We show be- 
low that in the experiment of Ref. || the " curtain effect" 
starts to work only at n > 50 mM. 

Two central parameters of our theory, A'o and e, can 
be calculated forra,siHiple model of DNA as an uniformly 
charged cylinderE.3~ll3. They can be also directly mea- 
sured in independent experiments. The energy e was 
measured in Ref. 

On the other hand, one can obtain these parameters 
from the experimental values of A'c and Nd- Below we use 
condensation and resolubilization for DNA with Spe'*+ 
to estimate e and Aq. We obtain e = 0.07 fc^T and 
A^o = 3.2 mM. TliG, first number reasonably agrees with 
experimental dataEll. 

Thus far we have talked about long DNA. For short 
DNA fragments, one should take into account the mix- 
ing entropy of DNA molecules which makes their solution 
more stable and makes the window between A^c and Nd 
smaller. We calculate the phase diagram for the conden- 
sation of DNA olygomers consisting of L/b bases {L is 
the length of olygomer, b — 1.7 A is length of double 
helix per phosphate). This phase diagram is shown in 
Fig. 2 on the plane {N, y) , where 
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C is the concentration of DNA molecules and Cmax is 
its maximum value (equal to the inverse volume of a 
DNA molecule). Below the solid curve, in segregation 
domain, condensed and dissolved phases coexist. Think- 
ing about very small concentrations of DNA, C ^ Cmax: 
we completely ignore the lower border of the segregation 
domain. Above the solid curve, all DNA molecules are 
in solution. A solution of long DNA molecules corre- 
sponds to a horizontal line y — 1. As concentration 
of DNA or its length decrease, the window between Ac 
and Nd shrinks. Condensation provides largest gain of 
free energy at A^ = Aq where DNA molecules are neu- 
tral. This is why the phase boundary curve peaks at 
A = Aq. The dotted line on Fig. 2 which obeys equation 
N = No devides the plane in two parts. At A^ < A^o 
the net charge of a DNA molecule is negative, while at 
A^ > A'o it is positive. This net charge by definition 
includes cations which arc bound to DNA with binding 
enery larger than ksT. Therefore, they move together 
with DNA in gel electrophoresis. This means that DNA 
should move in the direction opposite to the conventional 
one if A > A^o liZl. It would be interesting to verify this 
prediction experimentally. We propose also to combine 
experimental study of the condensation phase diagram 
with the gel electrophoretic measurement of the charge 
on the dissolved DNA molecules. 
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FIG. 2. Phase diagram of a dilute DNA solution in the 
plane of the bulk cation concentration A*' (in units of mM) 
and variable y defined by Eg. (|^). The solid curve was plot- 
ted using Eqsj-/^^ and ( |25| ) with fitting parameters obtained 
from the dataci for long DNA screened by Spe*^ . Above the 
curve, all DNA molecules are in solution {x — 1). Below the 
solid curve, the segregation domain is located. Here the con- 
densed phase of DNA appears and in regions far enough from 
the curve, it consumes most of DNA (a; <C 1). The dotted 
line corresponds to A'^ = A^o, where the net charge of DNA 
changes sign from negative to positive 



II. CRITERION OF EQUILIBRIUM BETWEEN 
CONDENSED AND DISSOLLVED PHASES OF 
LONG DNA. 

In this section we consider long DNA molecules in an 
aqueous solution containing Z-valent counterions of the 
bulk concentration N. DNA condensation and resolubi- 
lization take place when the chemical potential of DNA 
in its condensed and dissolved states are equal, i.e. 



Mc = fJ-d 



(2) 



For a long DNA we can neglect the entropy of DNA in 
both phases. We view the condensed phase as a bundle 
of parallel molecules which stick together due to correla- 
tion induced short range attraction to the nearest neig- 
bours. The chemical potential of DNA molecule in the 
condensed phase in this approximation is determined by 



fJ-c 



-eL/b, 



(3) 



where L/b is the number of negative charges in a DNA 
molecule. 

Each molecule of a large bundle is practically neutral. 
More exactly its charge is inversUt-, propotional to the 
number of molecules in the bundlellJ. This happens be- 
cause the total charge of a large bundle keeps counterions 
inside the bundle very effectively. 

On the other hand, in the dissolved state, each 
molecule aquires finite net charge density, rj*. This hap- 
pens because a finite fraction of its counterions moves to 
the Debye-Hiickel atmosphere in order to increase their 
entropy. One can view a single DNA molecule and its 



Debye-Hiickel atmosphere as a cylindrical capacitor with 
linear charge densities ij* (inside) and —r]* (outside). We 
recall that the net linear charge density of DNA, 77*, in- 
cludes the bare charge and the charge of cations residing 
at the very surface of DNA and attached to the surface 
with binding energy larger than ksT. 

In equilibrium the electrostatic potential on the sur- 
face of DNA, ipo, is determined by the ratio of the cation 
concentration in the bulk, N and the concentration near 
the surface of DNA, Nq: 



(4) 



This potential acts like a voltage difference applied to the 
capacitor. 

The free energy per DNA molecule or the chemical po- 
tential of the dissolved phase, can be written as 



1 *T 

/id = --^ifoV L. 



(5) 



To derive Eq. (^) one should add the energy of the electric 
field of the capacitor. 



(6) 



to the change of the entropy term of the free energy of 
cations when they move from the surface of DNA to the 
bulk of solution 



Lv* No 

— ksT ln- = -^or/L 



(7) 



(here Eq. (||) was used). Equivalent derivation can be 
found, for example, in Ref. |2^. Comparing Eq. (|^) and 
Eq. (^) we see that the change of the free energy of the 
capacitor is equal to the energy of an electric field with a 
minus sign. This is a realization of the general-theorem 
valid for any capacitor kept at constant voltageES. 

The surface potential, tp^, can be easily related to the 
net charge density 77*. Indeed, at distance r from its sur- 
face, a cylinder of radius a and linear charge rf creates 
a potential 



277* , I a + r 

In 

D 



(8) 



where D is the dielectric constant of water. This poten- 
tial vanishes beyond the Debye screening length. 



r, = {AnlBy^^^ {NZ^ + ZN + 2n) 



-1/2 



(9) 



where Is = /DksT is the Bjerrum length. Therefore, 
substituting r ^ and (p{rs) ~ into Eq. (^ one ob- 
tains 



277* fa + r, 
U \ a 



(10) 
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Substituting Eq. (|l^) into Eq. (|) we get 



(11) 



Using the condition (||) together with the expressions for 
the chemical potentials for both states, Eqs. (||) and ( |ll| ) 
we arrive at the final equation for the DNA reentrant 
condensation transitions 



knT 



D 



ln(l + r,/a). 



(12) 



To proceed further one has to know the net charge den- 
sity, 77*. Next section gives a review of current under- 
standing of this quantity. 



III. NET LINEAR CHARGE DENSITY OF 
SCREENED DNA 

Conventional understanding of nonlinear screening of 
a strongly charged cylinder is based on the Poisson- 
Boltzmann equation. Let us first consider a cylinder 
with a negative bare linear charge density, —rj, which 
is screened by Z-valent cations (n = 0). Assume that 
77 > Jjc, where rjc = c/Ib (j] =i-4-2 rjc for double he- 
lix DNA). Onsager and Manningll argued that such a 
cylinder is partially screened by counterions residing at 
the very surface, so that the net linear charge density 
of the cylinder, 77*, is equal to the negative universal 
value —"qc/Z. The rest of the cylinder charge is linearly 
screened at much larger distances according to the lin- 
ear Debye-Hiickel theory. The net charge, 77*, does not 
depend on the bulk concentration of cations, N and is 
shown (in units of — 7/c) by the dotted line on Fig. 3a. 
The Onsager-Manning picture of condensation was con- 
firmed by the solution of the Poisson-Boltzmann equa- 
tiorO. 

Let us now discuss a water solution, containing both 
a concentration, n, of monovalent salt and a concentra- 
tion, of Z : 1, salt. When N grows to some well 
defined concentration, N = N^, multivalent cations re- 
place monovalent ones at the swrfaee of DNA. According 
to conventional understandinyiri this replacement of 
the condensed monovalent ions changes ~Tj* /rjc from 1 
to 1/Z (see the dotted line on Fig. 3b). In logarithmic 
scale this transitioH-Joaks quite abrupt. 

In recent paperal3'tll the influence of strong correla- 
tions of multivalent cations at the surface of a macroion 
on its net charge, rf , was studied. The general expression 
for the net linear charge density, rj* , was derived: 



77 



Tj^ ln(Ao/7V) 
2Z \n{l + r,/a)' 



(13) 



This equation takes into account correlations with the 
help of the new boundary condition A^(7' = 0) = A^o for 
the Poisson-Boltzmann equation. To remind derivation 
of this result let us write down the Boltzmann formula 



N{r) = No exp 



Ze 



{ip{r) -ifio)] , 



(14) 



At r = Ts the concentration N{r) reaches its bulk value, 
A^. Substituting Eq. (0) into Eq. (^ we obtain 



f2Zrj\ fa + Ts 
N = Nq exp I In 



(15) 



Expressing 77* in terms of the bulk cation concentration 
N we arrive at Eq. (13). This equation was used to plot 
solid lines on Fig. 3. 
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FIG. 3. The dimensionless net linear charge density rj* /rjc 
as function of Z-valent cation concentration A'^ (in units of 
mM) at zero concentration n of monovalent salt (a) and at 
finite n (b). The solid curves are drawp according Eq. ( |l^ ) 
with parameters obtained from the datad for DNA screened by 
Spe''+. Dotted curves represent conventional understanding 
of results of Poisson-Bolzmann equation. A^o is the concen- 
tration of Z-valent ions at the surface of DNA and Nr is the 
concentration N at which Z-valent cations replace monova- 
lent ones. 

At 71 = (Fig. 3a) one can see0 that Manning's limit- 
ing value of 77* = —rjc/Z holds only for the unrealistically 
small A^ . At larger A^ the absolute value of the net charge 
decreases, 77* crosses zero and becomes positive. 

At finite ti, the screening radius, r^, strongly decreases 
and at small A^ the net charge density, —rj*, becomes 
larger than rjc/Z (see solid line on Fig. 3b). It was shown 
in Ref. |l^ that when A^ goes to Nr, the density 77* be- 
comes as large as 77c, so that at the replacement point the 
net density, 77*, matches its standard value for monova- 
lent ions. At larger A^ the role of the monovalent salt is 
smaller and 77* is similar to that of tt, = 0. The density, 
77*, crosses zero at A^ = A^o and becomes positive. In the 
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next section we will use Eq. ( |13| 
and resolubilization thresholds. 



to find the condensation 



IniV, = IniVn 



ln(l + ri/a) 



1/2 



(19) 



IV. REENTRANT CONDENSATION 
THRESHOLDS FOR LONG DNA MOLECULES. 
COMPARISON WITH EXPERIMENT 

Now we can substitute Eq. ( p^ into Eq. ([l^) and ob- 
tain an explicit equation for the threshold concentrations: 



1 



^ In' (Nq/N) 

kgT AZ^C ln(l + r,/a)' 



(16) 



Here ^ — Is/b is the conventional Manning's parameter 
(^ = 4.2 for DNA). Obviously, there exist two solutions, 
Nc and Nd, for the bulk concentration, N, which corre- 
sponds respectively to the condensation and dissolution 
transition points. 

Concentration is usually so small that is dom- 
inated by monovalent ions and = ri = {Snlsn)^^/^ . 
With the help of Eq. ([l6| ) we arrive at a simple equation 
for Nr-. 



e ^ ln'(iVo/jVc) 
fcsT 4^2^ ln(l + ri/a)' 



(17) 



To obtain an equation for we take into account that 
in the case of Spe''+ catioro, the condensed phase dis- 
solves at Nd ~ 150 mM. At such a large concentration, 
rs is determined by multivalent cations. However, if we 
use Eq. (^), rs turns out to be much smaller than the 
average distance between cations in the bulk of solution, 
which does not make physical sense. This means that the 
Debye-Hiickel theory does not work. A natural approx- 
imation in this case is to replace Ts by average distance 
to closest cation = (47rA"d/3)^^/'^. This gives for Nd 



1 In' [No/Nd) 



keT \Z\ ln(l + rrf/a)' 



(18) 



In practice, the concentrations N^. and Nd are known 
from experiments, so that important parameters e and 
A^o can be found with the help of Eqs. (^7|) and (18). 
For example, Nq can be found by eliminating e from 
eqs.(|l2, |l^). In the specific experimental situation of 
Ref. Wwhere long DNA molecules condensed with a 4- 
valent cation, Spe^+, N^ = 0.025 mM and Nd = 150 mM 
at the lowest concentration n = 10 mM of NaCl. Using 
these values we obtain from Eqs. (0), (|l|) Ao = 3.2 mM 
and e/ksT « 0.07. The last value favorably agrees 
with the energy of attraction between the CoHex cop. 
densed DNA obtained in the osmotic stress experimenttil 
e = 0.08 ksT. 

Knowing Nq and e, we can try to reproduce the exper- 
imental dependence of the condensation threshold N^. on 
concentration of NaCl, obtained in Ref. ^j. This can be 
done with the help of Eq. (0): 



The calculated function Nc{n) is shown on Fig. 4 to- 
gether with the experimental points from Ref. ^. It 
is clear that Eq. (|l^) closely reproduces experimental 
behavior till n < 50 mM. The later value is close to 
the concentration of monovalent cations which is needed 
to replace 4-valent cations at the surface of DNA for 
N « 0.01 mM. It seems that at n > 50 mM, conden- 
sation happens simultaneously with the replacement of 



monovalent cations by multivalent ones, i.e. at Nc 
as was assumed in Ref. E5. 
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FIG. 4. The condensation threshold as a function of 
monovalent salt concentration, The solid curve corresponds 
to Eq. (p^). Experimental dataa are shown by the large dots. 

We would like to emphasize that in the broad range of 
lower monovalent salt concentrations n < 50 mM when 
multivalent cations cover the DNA surface, condensation 
does not happen. In this range, the Coulomb repulsion 
of negative DNA molecules is strong enough to prevent 
condensation because, as shown on Fig. 3b, the absolute 
value of the net negative charge on the polymer, — ?7*, 
does not abrubtly decrease from rjc to rjc/Z at Nc Nr, 
but still stays larger than rjc/Z in the large interval of 

N > Nr. 

In contrast with Nc, the resolubilization threshold Nd 
practically does not depend on n. This happens because 
at A = Nd, screening is dominated by multivalent cations 
and agrees with our Eq. (p^. 

Let us now discuss another experimental studyQ in 
which DNA were condensed with trivalent Spd'^"'" cations. 
In this case, for rt = 4 mM thresholds are equal Nc ~ 
2 mM and Nd « 60 mM, which yields Aq w 11 mM 
and e/kBT « 0.02. This value of e/ksT is about three 
times smaller than for Spe*^. This is quite reasonable, 
if attraction has an electrostatic nature. 

The concentrations, Nq, we obtained abovjg^ are sub- 
stantially larger than the microscopic theoryli3 predicts 
for the point-like 3- and 4-valent cations. Actually, 
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cations Spd'^+ and Spe^"*" are quite long (« 15 and 20 A 
respectively) linear polymers. Their length is approx- 
imately equal to the average distance between cation 
centers in two-dimensional strongly correlated liquid of 
cations on the surface of DNA. Thus, to obtain a more 
reliable theoretical prediction of Nq, one has to study 
the thermodynamic properties of the strongly correlated 
(possibly nematic) liquid of these ions on an uniform neg- 
ative background by numerical methods. We will address 
this problem in our next paper. 

The application of our theory to experimental dat£ 
is not completely convincing, because of the following. 
Our theory assumes that the concentration of DNA is so 
small that the concentration of DNA phosphates Nph is 
smaller than N, so that the DNA charge can always be 
compensated by-Z-valent cations. Actually in experimen- 
tal conditionsO'cl the condensation threshold Nc happens 
to be close to Nph- It is possible, therefore, that Nc in 
this case is determined by the condition Nc ~ Nph rather 
than the repulsion of DNA helices. This means that the 
actual value of Nc can be lower than the measured one. 
Therefore, the true values of A'o should be somewhat 
smaller, and that of e should be somewhat larger than 
our estimates. We suggest a repetition of experiments at 
smaller concentrations of long DNA. 

Another way to improve the reliability of the extrac- 
tion of the value of Nq from experiment is to study the 
condensation of short DNA (see next section). 

Note that Ref. || also contains information on the sin- 
gle and triple stranded DNA helices, which condensed in 
the narrower and wider range of [Spe*+] than the double 
helices respectively. This tendency agrees with idea that 
correlations play the major role in this phenomenon. 



V. CONDENSATION OF SHORT DNA 
MOLECULES. 

In this Section we explicitly deal with the mixing en- 
tropy of DNA molecules in the dissolved state, and cal- 
culate the phase diagram for the DNA solution. We con- 
sider double helix DNA molecules of length L, with L/b 
bases each. 

The free energy per DNA molecule in solution with 
concentration C is 



1 Iri^iNo/N) 



F = -xkeT In 



Cx 



l|-(l-a;)^A, (20) 



where 1 — a; is the fraction of DNA molecules in the con- 
densate. The first term in Eq. ( pO| ) is the entropy of 
the dissolved DNA phase per molecule of DNA. Here 
Cmax ^ l/vra^L is the inverse volume of the DNA 
molecule. The second term in Eq. ( po[) is the average 
free energy per molecule in the condensed state. Here 
A is the difference between the energy of short range 
attraction e and the free energy bjid/L of the screening 
atmosphere per nucleotide (see Eq. (O)): 



keT ksT 4Z2^1n(l +rs/a)' 
Minimizing Eq. ( [20| ) with respect of a;, we find 

A L" 



(21) 



^max 

X = — — — exp 



C 



keT b 



(22) 



Then the boundary of two-phase domain x = \ corre- 
sponds to the condition 



In 



A L 



C 



kuT 6' 



(23) 



which can be explicitly written for N < Nq as 



bksT Cnax _ ^ ksT 1 \n\No/N) ^^^^ 



Le 



C 



and for N > Nq 
bksT , C„ 



Le 



■In 



C 



= 1 - 



e 4Z2^1n(l + ri/a)^ 



ksT 1 ln^(iVo/A^) 
^~4Z2^ ln(l + rw/a)' 



(25) 



where r^v = (47rA^/3)~^/'^. The last three equations spec- 
ify the conditions when the chemical potentials of the 
condensed and dissolved states of DNA are equal. The 
curve described by Eqs. ( p^ and (p5|), forms a bound- 
ary on the phase diagram, separating the region of the 
dissolved DNA (above it) from the region where the con- 
densed and dissolved phases of DNA coexist (below it). 
Within the separation domain the fraction of dissolved 
DNA is given by Eq. (||). 

The position and the shape of the phase boundary 
obviously depends on the values of the parameters Nq, 
e/ksT and 2Z'^£,. For double-helix DNA with Spe4+ we 
find 2^2^ = 134 and using the estimates e/fcsT = 0.07 
and Nq « 3.2 mM obtained above we can calculate the 
phase boundary for condensation of the DNA fragments 
of the arbitrary length. This boundary is shown on Fig. 2 
by the solid line. 

The dotted line N — Nq divides the phase diagram of 
Fig. 2 into two regions where the DNA molecules have 
opposite signs. Positive net charge of dissolved DNA can 
be measured in an electrophoresis experiment because 
cations included in this charge are bound to DNA with a 
binding energy larger than fc^T and therefore move to- 
gether with DNA molecule. At > A'o one should see 
that DNA molecules are positive both above and below 
phase boundary of Fig. 2. However, below the boundary, 
intensity of corresponding electrophoresis peak should 
decay rapidly with the distance from the boundary. This 
intensity is picked up by slowly moving bundles of DNA 
molecules, which at large enough N can be also positive. 
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The theory of DNA condensation and resolubihzation 
by multivalent cations presented above makes several 
novel, well-defined predictions which have not been con- 
firmed by experiment yet. The main results of the study 
are summarized in the phase diagram presented on Fig. 
2. There are only two physical parameters e and iVo on 
which the shape of the phase boundary depends. There- 
fore measuring just a few threshold concentrations of 
multivalent cations Nc and Nd for solution of DNA frag- 
ments of different length and/or concentrations should 
yield several independent determinations of these quan- 
tities, and at the same time provide the test for the self- 
consistency of our model. Experimental studies of DNA 
phase diargam with different multivalent cations would 
provide the values of the attractive energy for different 
ions. In the present theory the origin of parameters e 
and iVo was not specified and they were treated as phe- 
nomenological parameters. Comparison of experimen- 
tally determined e and iVo values for different counterions 
should yied information about the nature of the attrac- 
tion. It is worth noting here that if the attraction is of 
electrostatic correlation origin the value of e in principle 
includes attraction due to the correlations of multivalent 
cations as well as all of the repulsive non-Coulomb DNA- 
DNA interaction. However, the first interaction decays 
slower than the second one. Then the binding energy e is 
determined by correlations. In this^ca^e, the two quanti- 
ties e and Nq are not independentlljO namely In A^o oc e. 
For the model of the uniformly charged cyJAHder both 
parameters were calculated microscopicallyll30. 

We would like to emphasize that the concentration A'o 
found in this paper plays extremely important role in 
any phenomenon related to screening of DNA molecules 
by multivalent ionst£lllZI. In this paper, we try to at- 
tract attention to the fact that A'o plays major role in 
electrophoresis, because the net linear charge densite.p£ 
DNA rj* changes sign at = A'q. It was predictedllj'EZI 
that DNA should start moving in the opposite direction 
at A^ > A'q. It is not obvious that one can sec this phe- 
nomenon for long DNA. Indeed, in large interval of con- 
centrations Nq < N < Nd most of long DNA molecules 
are condensed in low mobilty bundles while concentra- 
tions A^ > Nd may be difficult for experiment because 
of large dissipation of heat. Therefore, we suggest doing 
electrophoresis of a solution of short DNA fragments. In 
this case all DNA molecules have unconventional sign of 
mobility at smaller than A^^ concentrations. The phase 
diagram shown in Fig. 2 predicts good conditions for 
such an experiment. It would be very interesting to ver- 
ify predicted correlations between the reentrant conden- 
sation and unconventional electrophoresis. 
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